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We derive a dimensionally-reduced limit theory for an n-dimensional 
^ ■ nonlinear elastic body that is slender along k dimensions. The start- 

CN ! ing point is to view an elastic body as an n-dimensional Riemannian 

manifold together with a not necessarily isometric W^'^ -immersion 
in n-dimensional Euclidean space. The equilibrium configuration is 
^ , the embedding that minimizes the average discrepancy between the 

^ I induced and intrinsic metrics. The dimensionally-reduced limit the- 

ory views the elastic body as a fc-dimensional Riemannian manifold 
along with an isometric W-^'^ -immersion in n-dimensional Euclidean 
space and linear data in the normal directions. The equilibrium con- 
figuration minimizes a functional depending on the average covari- 
ant derivatives of the linear data. The dimensionally-reduced limit 
is obtained using a T-convergence approach. The limit includes as 
particular cases plate, shell, and rod theories. It applies equally to 
"standard" elasticity and to "incompatible" elasticity, thus includ- 
ing as particular cases so-called non-Euclidean plate, shell, and rod 
theories. 
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1 Introduction 

The derivation of dimensionally-reduced elastic theories is a longstanding 
theme in material sciences, which goes back as far as to Euler, D. Bernoulli, 
Cauchy, and Kirchhoff | JKir50y , and in the last century, to name just a 
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few, to von Karman [ jvonl(J| ], E. and F. Cosserat, Love [ JLov27[ j, and Koiter 
[ |Koi66| j. In essence, the problem is the following: Elasticity theory models 



the static and dynamic properties of three-dimensional bodies. Unless 
simplifying assumptions can be made, these models are highly nonlinear 
and notoriously complex. In many cases, however, one is interested in 
elastic bodies that are slender across one or two dimensions. In such cases 
it is appealing to model the body as an object of lower dimension — either a 
surface or a curve, depending on the number of slender axes. The challenge 
is to derive reduced models for surface-like or curve-like bodies, departing 
from the full three-dimensional model. 



This is the context to which belong plate, shell, and rod theories [ frWK59 



Koi66[ j. Plates and shells are thin elastic sheets that are modeled as two- 
dimensional surfaces; plates are in a state of mechanical rest when they are 
planar, whereas shells are at rest in a non-planar configuration. Rods are 
thin and elongated bodies that are modeled as one-dimensional curves. 

Elastostatics can be formulated using a variational approach, where the 
equilibrium configuration of the body is a minimizer of an energy func- 
tional defined on the space of configurations. Until quite recently, the 
derivation of dimensionally-reduced models from three-dimensional elas- 
ticity had been based mostly on ansatzes for energy minimizers (e.g., the 
Kirchhoff-Love assumptions [ |Kir50| , ^ov2^]). Such an approach is non- 



rigorous, and in particular, different ansatzes lead to different reduced 
models (sometimes mutually inconsistent), a situation that has lead over 
the years to numerous controversies. This situation has changed drasti- 
cally over the last decade with the development of new analytical methods, 
based on T-convergence, which have lead to rigorous derivations of plate 
[ |I-'JM02[ J, shell I 1EJMM03I J, and rod | |MM03[ | theories. 



An underlying assumption of classical three-dimensional elastic theories 
is the existence of a configuration in which the body is free of any inter- 
nal stresses; this reference configuration is unique modulo rigid transforma- 
tions. Deviations from the reference configuration involve an elastic energy 
"cost", and only occur in response to external forces or to the imposition 
of constraints, such as boundary conditions that are incompatible with the 
reference configuration. In the absence of such forces or constraints, the 
reference configuration is the equilibrium configuration. 

There are many systems however, in which a stress-free configuration 
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does not exist; such bodies are said to be residually-stressed. The distinctive 
feature of a residually-stressed body is that its constituents change their 
shape if the body is dissected into parts (thus releasing the residual stress). 
The study and modeling of residually-stressed bodies has its own history. 



starting with the pioneering work of Wang [1Wan66i , |Wan670 and Kroner 



[ |Kro81| j in the context of defects and dislocations. In more recent years 



there has been a renewed interest in residually-stressed elastic bodies in 
the context of pattern formation in plants (see e.g. [ |VGU8| , |LMU9| , |AESKll| j) 



and in synthetic materials, such as thermo-responsive gels QKESO/I ]. 

A residually-stressed body (made of an amorphous material) can be mod- 
eled as a three-dimensional Riemannian manifold M; the intrinsic property 
of the material is (local) distances between neighboring material elements. 
Thus, the concept of a reference configuration is replaced by that of a refer- 
ence metric q\^^^^. Aconfigurationof this body is a mapping/ : M — > 
from the Riemannian manifold into the ambient three-dimensional Eu- 
clidean space. The mapping / induces on M a metric g = f*^, where e is 
the canonical Euclidean metric in R". The local energy density associated 
with the mapping / is assumed to only depend on the metric discrepancy 
between g and g (this metric discrepancy is called a strain is the elastic 
context). 

If M is simply-connected and the reference metric g is flat, namely, has a 
vanishing curvature tensor, then there exists an isometric immersion / of 
(M, g) into (R", e), in which case g = Q, hence / is a stress-free equilibrium 
configuration. If g is not flat, there does not exist a configuration for which 
g = Q, and the energy minimizing configuration will necessarily carry non- 
zero elastic energy, i.e., it will be residually-stressed. In the elastic context, 
the metric g is said to be incompatible, and the elastic model is known as 
incompatible elasticity. 

Organisms such as plant organs, and the thermo-responsive gels studied in 
[ |KES07| ] are surface-like residually stressed materials. The interest in such 
systems has naturally lead to the development of dimensionally-reduced 
models in incompatible materials. A reduced theory of non-Euclidean plates 
was developed in [ [EbK(J9bP based on a non-rigorous ansatz, followed in 
[ |ESK09c| ] by a reduced theory of non-Euclidean shells (we will elaborate 
on the distinction between the two cases in Section H). Lewicka and 
Pakzad [ p^PlUI ] rigorously derived the limit of non-Euclidean plates us- 
ing a T-convergence approach, thus generalizing the result of Friesecke et 
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al. [ |bJMU2| ] (the latter being a particular case for a flat g). 



In this paper we present an analysis that generalizes in one fell swoop the 
derivation of (Euclidean) plate | jJ-'JMU20 , she ll Ji-'JMMOl , and rod [ |MMU3| ] 
theories, and (non-Euclidean) plate theory OLPIOH , and also provides, as 
particular cases, a rigorous derivation of non-Euclidean shell and rod theo- 
ries, for which no current analyses exist. Specifically, we consider a family 
of elastic problems in which the domains Q;, c M are a one-parameter 
family of n-dimensional Riemannian manifolds that are shrinking into an 
(n - ^:)-dimensional submanifold S as the thickness parameter h tends to 
zero. For each such domain we consider immersions into an n-dimensional 
Euclidean space (in the physical context n = 3 and k = 1,2) and associate 
with each such immersion an energy : W^'^(Q;j;R") — > R. Under suit- 
able assumptions, we prove that any family of (possibly approximate) 
minimizers fh s W^'^{D.h', IR") of converges in a sense we define (modulo 
subsequences) as /z — > to a mapping / 6 W^''^(S;R"), which minimizes a 
limiting energy functional £iim : W^'^(§;R") — > R. Thus, the elastic problem 
associated with the immersion of the (n - fc)-dimensional manifold S into 
R" is the dimensionally-reduced model in the h ^ limit. 

As stated above, this general result embodies the existing theories for 
plates, shells, and rods, the distinction between the first two cases being 
solely a property of the reference metric g, whereas the distinction between 
plates/shells and rods is the codimension k. Indeed, an approach based 
on a reference metric rather than a reference configuration makes the dis- 
tinction between plates and shells almost unnoticed. Moreover, under this 
viewpoint, there is nothing special about non-Euclidean plates, shells, or 
rods, except for the fact that the limiting energy functional £iim may not 
assume a state of zero energy. The remarkable fact is that the formulation 
of the problem within the framework of Riemannian geometry leads to a 
limiting model that has the exact same form in all instances. 

While our result is very general in that it covers a variety of limits of elastic 
problems that were previously treated separately, we have deliberately re- 
stricted our attention to a specific energy density, which can be considered 
as a generalization of Hooke's law for isotropic materials, and to slender 
bodies with a symmetric cross section. An extension of the present analysis 
to remove these restrictions seems to be straightforward. 
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2 Problem statement and main result 



Let M be an n-dimensional smooth oriented manifold; let § c M be an m- 
dimensional smooth oriented submanifold; let k denote the codimension 
of S so that m + k = n. We endow M with a smooth Riemannian metric g. 
The submanifold (S, g|s) is bounded, either closed, or having a Lipschitz 
continuous boundary. In either case, the Riemann curvature tensor is 
uniformly bounded in M. In the context of incompatible elasticity, the 
manifold (M, g) models an elastic body with internal distances prescribed 
by the reference metric g. 

Let h 6 (0, ho) be a continuous parameter, and let 

Q,, = {p G M : dist(p,S) <h}cM 

be a family of tubular neighborhoods of S that inherit the reference metric 
g. The boundedness of the curvature implies that for small enough h, the 
exponential map, 

exp:{(p,e)GKS: 1^1 < ^ Q, 

is a diffeomorphism, where 3\fS = TM|g is the orthogonal complement of 
TM\l ^ rS in TMIs. Thus, we identify 

TM|s = TS e XS. 

As standard, we denote by n : !NS — > S the projection from the vector 
bundle NS, or its restriction to Q;,, onto its base. Let E — > S and f — > 3Nf§ 
be vector bundles, and let O : n*E — > f . Let £, G !N§ and let rj G £71(5) • 
The fiber (7i*E)^ is canonically identified with the fiber So, we can 

unambiguously apply O at ^ to r/. Denote the result by 0^(7]). 

We next introduce some definitions and notations. We define the projection 
operators, 

P" : TMIs ^ TS and : TM\§ XS, 

and the corresponding inclusions 

i" : rS ^ TM\s and : K§ TMIs- 
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Let V be the Levi-Civita connection on TMIg. Then, the induced connec- 
tions on TS and !NS are 



V" = P" o V o i" and = o V o l^. 

When it does not cause confusion, we use V without any decorations to 
denote the induced connections as well. 

Let L : 7i*!NS TJ^S denote the canonical identification of the vector 
bundle JsfS with its own vertical tangent space. Specifically, for ^ G KS 
and 7] G (7i*!NS)^, there is a canonical identification of i] with an element of 
(N§)„(^). We then define a curve y : I ^ ?\f§, 

y{t) = ^ + rit, 

and identify L^irj) = y(0). 

As stated above, we identify Q;, with a subset of 3\fS. Consider now the 
tangent space TIN'S. Define the isomorphism 

n : nT§ e nJ<§> TKS, 

as follows. Let C denote the zero section of !N"S. Define Tl to be the unique 
map such that 

nis = dc®4, 

and for each ^ G !NS we have 

v,,(£)n = 0. 

That is, n is given by radial parallel transport along the fibers of XS. 
A notable property of n is that it preserves the metric, namely 

Q{n^u,n^v) = q{u,v), g XS, t/,t; g Tn(£)M. 

For every h G (0,/zo) we consider mappings G W^'^(Q;,;R"), and assign 
to every such mapping an energy 6^ [/;,]. In (hyper-)elasticity the energy 
is assumed to be a volume integral of a non-negative energy density Wh 
[ [Tru55[ J. In the absence of external constraints and forces, W;, only depends 



on the local value of the derivative of the mapping dfh, and only vanishes if 
is a local orientation-preserving isometry, namely, if d/;, G SO(n), where 

SO{n) = {q : TQ;, — > R" : q*e = g, q is orientation-preserving} . 
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The mappings G W^'^(Q;,;R") are only of interest modulo rigid trans- 
formations, hence we may assume (for the sake of a later compactness 
argument) that 

r fHdvol, = 0. (2.1) 

In the present work we consider a specific energy functional that postulates 
that the material is isotropic and that the energy density scales quadrati- 
cally with the distance of dfh from SO{n). Such an energy density can be 
viewed as a continuum variant of Hooke's law for linear springs (i.e., the 
energy density is quadratic in the local strain). Specifically, we define 

^h[fk] = ^£ dist'(d/,, SO(n)) rfvolg, (2.2) 

where J denotes integration divided by the volume of the domain, and the 
additional pref actor is discussed next. 

For every fixed h G {0,ho) the energy functional ( |2.2| ) defines an elastic 
problem: find the mapping G W^'^(Q;,;]R") that minimizes £;,. It is not 
known a priori that such minimizers do exist, but one can always consider 
a family /;, of approximate minimizers, defined by the condition 

lim (£„[/„] -£;) = 0, 

where 

£;=inf {£,[/]: / G Wi'2(Q,;]R")) . 

As — > the family of n-dimensional domains Q;, shrinks to the m- 
dimensional submanifold S, and hence the volume integral of the energy 
density in the (exact or approximate) equilibrium configuration is expected 
to tend to zero. Since we are interested in the — > limit of this family 
of elastic problems, we have first rescaled the energy by dividing it by the 
volume of the domain. Furthermore, motivated by the physical setting 
in which S is a either a one- or two-dimensional submanifold of a three- 
dimensional manifold, one can expect Q;, to be "almost" W^'^-isometrically 
immersible in R", in the sense that even the energy per unit volume tends 
to zero as ^ — > 0. This amounts to the submanifold (S; g|s) being isomet- 
rically immersible in R" more regularly than W^'^. As will be shown, if 
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(S;g|s) is W^'^-isometrically immersible into R", then the energy per unit 
volume is OQi^) as ^ — > 0, which is why we divided the energy per unit 
volume in ( p..2\) by the additional 1 /h^ factor. 

Our assumptions about the immersibility of (S, g|s) in R" are encapsulated 
in the so-called^inzte bending assumption: 

There exists a sequence of mappings G W^'^(Q;,; R") such that 

^h[fh] = 0(1). (2.3) 

We now state our main results: Let 

eT = {(f, q^) : f G W^'\§; R"), G W^'\§; XS* (8 R")) • 

We say that a sequence of maps fh : O/i — > R" reduced-converges to an 
element (f, q^) G if 

lim f \fh-Fo 7zp rfvolg = 0, (2.4) 

and 

lim f \dfh on- 71* (if e q^)p dvolg = 0. (2.5) 

For the physically-oriented reader, condition ( |Z^ ) states that the confor- 
mations fh of the shrinking domains Q;, converge in the mean-square to 
a conformation f of the submanifold S. Condition ( |275| ) states that the 
tangential component of rf/;, consistently converges to the derivative of F, 
whereas the normal component of dfh converges to a limit q-"-. 

Here and throughout this paper, we denote by n*dF the section of 7i*r*S(8)R" 
obtained by pulling-back dF considered as a section of T* S(8)R" . This should 
not be confused with the closely related pull-back of dF considered as a 
1-form on § involving composition with dn, which we denote by 

n*dF = ndF o dn = dn o ndF. 

The first step in our analysis is to show that any family of (possibly ap- 
proximate) minimizers of reduced-converges, modulo subsequences, 
to an element (f, q-^) of the space ^ . 



9 



We introduce a functional £iim : ^ — > R U {00} , defined as follows. For 
(f, q^) G ^, define q G ^^'^(S, rM|s ® M") by 

q = rff ® q-^, 

and let G r(S; R" (8 TS) and G r(S; R" ® KS) be defined by 
pj = pll o q-i and P^ = P^ o q-^. 

Define 



£iim[f, q^] = 



^ jfl (2|P; o Vq^ - II)|2 + |P^ o Vq^p) rfvolg,, q G SO(n) a.e. 
1 00 otherwise, 

(2.6) 

where g|s is the induced metric on S and k is the volume of the k - 1 
dimensional unit sphere divided hy k + 2. 

We prove the following: 

1. The lozver-semicontinuity property, 

liminf £;,[/;,] > £iim[f,q-^]. 

2. For every (O, p) G there exists a family of mappings (/);, G W^'^(Q;,;R") 
(a so-called recovery sequence), such that (ph reduced-converges to 
(O, p) and 

limE,h[(ph] = £iim[0,p]. 

It is easy then to show (see e.g., dal Maso [ |dal93[ j on F-convergence) that 
the (possibly partial) limit (f, q-^) of the sequence /;, of (possibly approxi- 
mate) minimizers is a (true!) minimizer of the limiting functional £iiin/ and 
moreover that 

£iim[F, q"^] = lim £;,[/;,]. 

The practical implication of this result is the following: whenever faced 
with the need to find a minimizer : Q;, — > R" of £;, for sufficiently small 
h, one can rather look for a minimizer (f, q-"-) of £iim/ which approximates 
fh in the sense of ( |2.4D and ( [2.51 ). In most cases, the latter task turns out to 
be easier. 
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3 Preliminaries 



In this section, we collect a number of definitions, facts, and basic lemmas 
used throughout the paper. 



3.1 Integration 

Let M be a manifold with boundary and let N be a manifold without 
boundary. Let / : M — > N be smooth. Let E — > N be a vector bundle. 
Denote by A'(N, E) differential forms of degree I on N with coefficients in E. 
Denote by /* : A*(N, E) A*{M, fE) the pull-back of differential forms by 
/. If / and f\dM are proper submersions of relative dimensions k and k - 1 
respectively, denote by /* : A*{M,f*E) — > A*{N,E)[-k] the push-forward 
operator or integration over the fiber of /. See Bott and Tu [ pT82| ] for a 
discussion of integration over the fiber in the case when M is the total space 
of a vector bundle, which is what we will use. The push-forward operator 
fi, satisfies the following properties: 

1. Let N be the point so that E, fE, are trivial bundles. If a G A'(M, fE), 
then 

\ La / = dimM, 
fa = \-'M 

1 otherwise. 

2. If |6 G A*{M, fE) and a G A*{N, E), then 

A|S) = a A/,jS. (3.1) 

This is a generalization of the linearity of integration to the fibered 
context. 



3. Let 




be a commutative diagram of smooth maps, where / is a proper 
submersion, P is the fiber product M Xn Q, and g, h are the canonical 
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projections. Then, h is a proper submersion, and if a e A*{M,f*E), 
then 

Kg*a = k*fi,a. 

This is a generalization of the classical change of variables formula. 
It is easy to see that properties uniquely characterize Moreover, 

4. Let 

P ^ M N, 
where g and / are proper submersions. Then 

/* ° = (/ ° gh- (3.2) 
This is a generalization of Fubini's theorem. 

5. Let V denote a connection on £ as well as the associated pull-back 
connection on f*E. Let co G A'{M,f*E). The following generalization 
of Stokes theorem holds: 

W{fM = /*(Va;) + {-lymbuUco). 

The following special cases will be particularly useful. It follows from ( [3.1[ ) 
that if « = F is a zero-form on S and jS = dvo\, then 

7I*(f o 71 rfVOlg) = f 7T*rfVOlg. (3.3) 

It follows from p.2| ) and item |1| above that if 

Qh S point, 
then for every differential form (5 on !NS, 

r (^) = f fi, (3.4) 

Js Jq/, 
12 



3.2 The tangent bundle of a tubular neighborhood 



In Section ^ we defined an isomorphism IT : 7i*TS © 7t*!NS — > TKS using 
the connection on M and the identification of an open subset of 3\fS with 
a tubular neighborhood of S in M by the exponential map. In the present 
section, we will construct another isomorphism 

using only the induced connection on XS. Lemma below estimates the 
discrepancy between a ® l and Yl. Thus, we may take advantage of the 
linearity of ct i to expedite calculations. 

Let L : n*J^§ ^ TJ^S denote the canonical identification of the vector 
bundle with its own vertical tangent space, as explained in Section ^ 
The differential dn : T!NS — > 7z*rS is also defined canonically. Clearly, 

dnoL = 0, (3.5) 

which implies that 

is a short exact sequence. 

To fully determine an isomorphism TIN'S = n*T§ © n*J^§ we need a map 

o : nT§ TKS, 

such that 

dn o o = ld. (3.6) 

Here, we use the induced connection on !N"S. Define o to be the unique map 
such that for any curve « : I — > S and any parallel normal field <^ : J — > 
along a, we have 

Thus, we have constructed an isomorphism, 

a®i: nT§ © nJ<§, TKS. 

Let A denote the tautological section of n*J^S. That is, for ^ G 3\fS, 

= 5 G (7I*XS)£. 

Let II : 3\fS (g) TS — > TS denote the second fundamental form. Then: 
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Lemma 3.1 We have 

o®L-n = oonllo{A® 7i*P'l) + 0(A®2). 

Proof: We start by writing 

c7ei-n = (cj-no 7z*i") e (t - n o ni-^). 

First, we show that for ^ G 

l(i - n o 7zV)d = o(|5p). 



(3.7) 



Indeed, let p G S and C/ s XSp be arbitrary. Consider the sections l{Q, t(r]) 
of 7i*!NS|n-i(p). Thinking of p as a point in KS, we claim that 



V,(c)i(r/)lp = 0. 
In fact, by symmetry of the connection. 



(3.8) 



_ D(d_ 
ds \ dt 

= V,o,)i(C)lp 



{tC + S7]) 



s=0) 



t=OI 



t=0 



s=0 



(3.9) 



Since the identification between D.^ and a neighborhood of S c !NS is via 
the exponential map, we have 



V.(Qi(C)lp = 



(3.10) 



for arbitrary C- Equation ( [3. 81 ) follows by the polarization identity from 
equations i^) and ( p^ . 

An immediate consequence of ( pi8| ) is that 



So, by definition of IT, 



l(V.(c)K/]))sd = 0(|s^|). 



|v,(,c)(i - no Osd = 0(1351)- 
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Since 

Ho i-^ls = 4, 
equation p77| ) follows by integrating radially. 
It remains to show that 

(T - n O 7T*f" = C7 O 71* II o(A 7I*P'I) + 0(A^2). 



(3.11) 



Let p G S, let 7] G TpS and ^ G KSp. Let « : I ^ S with d(0) = rj. Let j6 be a 
parallel normal field along a with (}{0) = ^. For s G R, SjS is also a parallel 
normal field along a. So, by definition of o, 



t=0 



Applying the covariant derivative V to both sides, we calculate how o 
varies radially along ^ : 





ds 




J 


s=0 


(3.12) 




-1 

dt \ 




s=o) 


t=0 





= i"oII(^,??), 

where the last equality follows from the definition of the second funda- 
mental form. On the other hand, by Leibniz's product rule. 



£c7,,(n(se,r])) 



s=0 



= i"oII(5,T]). 



(3.13) 



Combining equations ( |3.12|) and ( |3.13| ) we obtain 



£(7s4(t])Is=o = £c7s4- (II(s<^, r/)) 



s=0 



So, considering p as a point of JnTS, by definition of IT, 

V,p(,c) (cT - n O 71* - CJ O 71* II o(A ® 7T*P")) = 0. 

Since 

(c7 - n O 7Z*i")|^ = = C7 O 71* II o(A 7Z*P")|s, 

equation ( |3.11| ) follows by integrating radially. 
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3.3 The metric on a tubular neighborhood 

CoroUary 3.1 

g o (cT e if^ - 7i*g = 0(A). 

Proof: Since parallel transport preserves g, we have 

g o ^ 7i*g. 

So, the corollary follows from Lemma ■ 

Let g denote the unique metric on the total space of such that o ® lis 
an isometry. It's easy to see that g|s = g|s- The following Corollary follows 
immediately from the previous. 

CoroUary 3.2 

0-3 = 0(A). 

In the rest of this paper we will write g instead of 7i*g when it does not 
cause confusion. 

3.4 The volume form on a tubular neighborhood 

Let £,F — > M, be vector bundles and let ^ : E — > F be a morphism of 
vector bundles. Denote by Kh : A^E — > A"F the associated vector bundle 
morphism between the a^^ exterior powers of E and F. 

Taking the determinant of n, we have an isomorphism 
Write 

P = (c7 e 0"^* ° (h*p")* : 7z*T*s ^ rm, 

= (cj e 0"^* ° {nV^y : nm* TJ^S. 

Note that 

(T*op = Id, ro0 = ld, 0*0 = 0, L*op = 0. (3.14) 
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Moreover, equations ( [3.14[ ) uniquely characterize p and 6. Taking exterior 
powers, we have 

A'p : A'nTS ATJ^S, A'd : A'nm* AiVJ^S. 

Moreover, 

^ A'pAAW = A\a®L)-^* : ^ A'nTS^A'nJ^S* ^ A'n{rS®y<S*) A'TJ^S. 

i+j=l i+i=l 

(3.15) 

Let fj be the unit norm section of A^T*§ belonging to the orientation class, 
i.e. fj = dvolgig. Let d) be the unit norm section of A'^'NS* belonging to the 
orientation class determined by the orientations of M and S. Define 

7] = A"'p o 71*7], CO = A^d o 710). 

In particular, 7] G A"'(^S) and co G A''^(IN"S). It is immediate from the 
definition that 

rj Aco = rfvolg. (3.16) 



CoroCCary 3.3 



T] Aco - dvolci = 0(A). 



Proof: The Corollary is an immediate consequence of equation ( [3.16D and 
Corollary ■ 

Lemma 3.2 Let a e A\§). Then 

n*a = a'p o Tia. 



Proof: By equations ( p.5D and (p.6D, we have 

p = dn*. 

So, the lemma follows from the definition of n'' 



Lemma 3.3 



Tj = 7l*rfuoZg|g. 
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Proof: Combine Lemma for I = m with the fact that 

A'"p o 7I*dvolg|g = A"'p o 71*7] = 7]. 

m 

Lemma 3.4 We have 

- -7^ = Oih), \§\v,h'^ - \Qn\ = OQi'-'), (3.17) 

|LJ;,| \b\ 

where is the volume of the k dimensional unit hall. 

Proof: By Corollary and Lemma we have 

7i*.((i')dvolg|g - 7i*.dvolg = 7i*(a; A n*dvo\^\^) - 7i*dvolg = 0{h^^^). 
Note that 7T*(aj) is the constant Vkh^. Integrating over S, we have 

\$\v,h^-m = o{h^'-^). 

The lemma follows. ■ 

Lemma 3.5 We have cuI^q,, = 0. 

Proof: Let S;,?\f§ c denote the radius h sphere bundle inside of Ji$. 
Since we have identified Q.h with a subset of via the exponential map, 
it follows that (90^ = S;,IN"S. Let A""- c 7i*J\fS denote the rank - 1 subbundle 
that is the orthogonal complement of A in NS. By definition of l, 

Moreover, since the connection V used to define o is metric. 

Counting dimensions, we conclude that 

is an isomorphism. So, it suffices to show that 

Indeed, by the definition of co and the third of equations p.l4| ). 

But, = because A""- has rank k- 1. The lemma follows. ■ 

18 



3.5 Rescaling a tubular neighborhood 

Define [i,, : Q;,^ Qn^h by 

Clearly n o ^z;, = n. So, there is a canonical bundle map fih : 7i*TM|s — > 
7T*rM|s covering /j;,. By abuse of notation, we use fi], to denote its own 
restriction to the summands n*T§ and 7t*!N"S of 7i*TM|s. 

Lemma 3.6 We have 



dy.h o L = h{Lo fijj) 
d[,Lh o a = a o fijj. 



(3.18) 
(3.19) 



Proof: Let^ G Qfog,letC s !NSn(^) and let x £ Tn(£)S. By definition i^(C) = y(0), 
where y(0 = + 5- So, 



di^h o ii(C) = ^fi;, o y(0 



Furthermore, 



f ;i ° y(0 = ^^C + h^- 



On the other hand, Lh^{Q = 6(0) where 5(0 = iC + /J^. Equation ( P^ 
follows from the fact that 



= h6{0). 



Similarly, by definition o^{x) = v(0), where v{t) is a parallel normal field 
along a path a{t) in S with d{0) = x and v(0) = ^. So, 



dl^h O C7^r(x) = —[^h O V(0 



f=0 



Moreover, /j;, o v(i) = /zv(f). On the other hand, hv{t) is a parallel normal 
field along a{t) with /zv(0) = h^. So, 



and equation ( [3.19D follows. 
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Corollary 3.4 We have 



L o d^il = h{fil o r), a o di.ll = fil o o\ 
d^il oe = h{eo jll), d^il op = pofil. 



(3.20) 
(3.21) 



Proof: Equations ( [3.20D are the duals of equations ( |3.18| ) and ( |3.19|) . Equa- 
tions ( |3.21| ) follow from equations p.2(J| ) and the fact that equations ( [3.14p 
characterize p and 0. ■ 



Corollary 3.5 ]Ne have 



= h^(V, [1*7] = 7]. 



Proof: The corollary follows from the definition of co and rj along with 
equations ( |3.21| ). ■ 



Lemma 3.7 Let f G \}(0.h^]^. Then 



/ / 



dvoL 



1 + 0{h) 

Vk\§>\ 



I (/ ° A 7]. 



Proof: Using Corollary pl3| . Lemma and Corollary p3| , we calculate 

f . , 1 + 0(h) r , „ 



i + o()i) r ,^ , 



Let 

4 : A'7z*rS ^ A'"-'7z*rS, : Nnm* P^'m^* , 

denote the Hodge star operators induced by the metric g. Let 

: A'rKS ^ A"-'rKs 
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denote the Hodge star operator induced by the metric g. Then, 



5' 



Y^{po.'^oo*)^{eoJ^on. (3.22) 

/+;=/ 



Denote by ft** the pull-back Hodge star operator, i.e. 

So, jU*? is the Hodge star operator of the metric g. 
Lemma 3.8 

K* = *- 

Proof: The lemma follows from equations ( ^^221) , and ( ^2^ . ■ 

3.6 Weak convergence on shrinking tubular neighborhoods 

definition 3.1 A sequence ofL^ differential forms an e A^^C^Mo) weakly con- 
verges to zero if for all O e A"~X^ho) have 



f (K 



fls /z — > 0. A sequence of sections a^ G L^(Qw,(,, 7z*A'r*M|s) weakly converges 



to zero if the corresponding sequence = A'(n*) ^ o g A' (O/jfeg) weakly 



converges to zero. 

Lemma 3.9 A sequence ah G L'^{Q.hho, 7t*A'T*M|s) weakly converges to zero if and 
only if 



Q,, 



Q{ah, ^h)dvol^ 

for all sequences fih of the form 

Ph = y-p{n^ o 

where p g L^(S; AT*S ® A'KS* (8) K§*®P) with i + j = n-l and p is arbitrary. 
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Proof: Let 

In particular, is a family of (n - l)-forms on dhgh- It is clear from the 
definition of A that A o y.^ = hA. So, by equations ( [3.21D we have 

The right-hand side of the preceding equation is clearly independent of h. 
So, we write 

Let ah = A^(n*)"^ o By Lemma PTT] , Corollary pTT| , equation ( p.lSI ), Corol- 
lary ^]2| and Lemma ^771 

f g(fl;7, jS;,)rfvol3 = I g(A^(c7 ® 0* ° och, j6;,)rfvolg + 0(/z) 

= f g(a;„ (A'p ® A'd) o |S,)rfvolg + 0{h) 

Jq,, 



(l + 0(^))f §(a,,|S,)rfvolg + 0(7i) 
1 + 0{h) 



Q,, 



V,c|S| 



(fi*g)(iU;>„0^)a; A7] + 0(h). 



Observe that is an arbitrary form on 0;^^ that is polynomial along the 
fibers of n. The lemma follows since polynomials are dense in L^. ■ 



4 Rigidity 

The compactness property, whereby any sequence of approximate mini- 
mizers of Eh reduced-converges, is based on a rigidity theorem that can be 
viewed as a quantitative version of Liouville's theorem. A rigidity theorem 
for mappings R" — > R" was proved by Friesecke et al. [ [bJMU2| ], paving the 
way to their derivation of plates, shell, and rod theories. In this section, we 
present a generalization of the rigidity theorem of [ pJM02| l that applies to 



our Riemannian setting. Like Lewicka and Pakzad in HLPlUl j, we base our 



proof on the theorem in Euclidean space. The notable difference between 
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Figure 1: 



our formulation of the rigidity theorem and the ones in [ |FJMO20 and [^PTO|] 
is that in the Riemannian setting one has to adapt the notion of a spatially 
constant matrix. Another difference between the approach here and the 
above mentioned references is the use of a smoothing convolution operator 
rather than a partition of unity. 



We introduce some more notations. Consider the commutative diagram 
in Figure ||. Recall that n : — ^ S denotes the canonical projection. 
Moreover, let CD : TS — > § denote the canonical projection, and let e : TS — > S 
denote the exponential map. The other maps in the diagram are canonical 
projections of fiber products. 



Below, we will use repeatedly the following identities that follow from 
the commutativity of the diagram, and the properties of the push-forward 
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operators: 

( e* = I CO* (4.1a) 

Js Js 

7i*e*=e*7i* (4.1b) 

O^n^e* = e cd*ti* (4.1c) 

e*7T* = (4. Id) 

e*cD* = cD*e* (4.1e) 

^*(D* = ao* (4.1f) 

^^a = (D*(Di,. (4.1g) 

For every h G (0,/zo) we define the indicator function /;, : TS — > R: 

hip, t]) = l|;j|<;,. 

We then define a family of non-negative test functions Kh : TS — > R whose 
support is compactly embedded in the support of the and satisfying 

CD* (X;, 7T*e*rfvolg) = CD* (Kh e 7i*rfvolg) = 1, (4.2) 

where the first equality follows from (|4.1bD. We further choose such that 



C^h^K,, < h < C2h"K2„, (4.3) 

and 

Ch"^'m\ < h. (4.4) 

Next, we introduce the function O : TS — > TS defined by 

e o O = CD and cD o O = e . 

Clearly, O is a diffeomorphism, and 

O o O = Id (4.5a) 

o O = /„. (4.5b) 

Changing variables, we may also express the normalization of Kh as fol- 
lows: 

e* {{Kh o O) cO*7z*rfvolg) = 1. (4.6) 
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We then define the mapping H : e* TS — > cdTS: 

H(p, 7], C) = (p, T], t), 

where t g TpS is the unique vector satisfying: 

expp(T) = exp^^p^(,^)(C). 

We note that 

j6 o H = CD 

e oa; o H = e oe" 

that is, the following diagram commutes: 

e* rS TS — ^-^ S 



(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 



We augment this diagram by adding the maps H, a and n as follows: 



e e*KS 



a* e* y^§ 



CD* rs 



hence 







_ VI 




a 






n o ^ = ^ o n 
e o a o S - e oe 



(4.8a) 
(4.8b) 
(4.8c) 



Let ^ : (D*TM\s — > e* TMIs be the isomorphism given by parallel trans- 
port along geodesic rays; we can view ^(p,,,) as a mapping from TpMIs to 

Texp (j])^ls- 
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Lemma 4.1 



(4.9) 



Proof: Holonomy around a loop is the integral of curvature over a spanning 
surface. ■ 

We these preliminaries, we state a local rigidity theorem: 

■Theorem 4.1 There exists a constant C > such that for every h G (0,/zo) and 
every fij e W^'^(Q;,;R") there exists a section p;, 6 L^(S; r*M|s <8) R"), suc/z f/zfli 



0. 



71* |[e*d/;, o e*n o 7i*M^ - 7i*cD*p;,|^ e*dyoZgj| 
< C {cD* [/;,7T*e* (dist^(rf/;j, SO(n)) rfyoZg^j + h^ cD* ^//j 7i*e*rfyo/gj| . 

Proof: Fix p G S. Let B/jTpS c TpS denote the ball of radius h, and let 

For sufficiently small h, we identify Upj, with the open subset of Euclidean 
space 



Vh = \{xi,...,x„)eK" 



as follows. Let rf be a basis of TpS and let E,^ be a frame of e* 15,^7^5 such 
that V^^ vanishes at G TpS for all The map u^^^ '■ Upji given by 

^ ] 

Up,h{x) = Xj+mE,' Xilf 

j=l \ !=1 / 

is a diffeomorphism and, therefore, defines coordinates on Up^h- We claim 
that with respect to the coordinates x„ the metric has the form 

9,7 = 6ij + 0{h^). (4.10) 

We return to the proof of ( |4.1U| ) below. Using this system of coordinates, we 
view maps Up^h R" as maps R" — > R". By the rigidity theorem proved in 
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[ |bJMU2[ J, there exists a constant C > such that for every / G W^'^{Upy, W) 
there exists an n x n matrix Q G SO{n), such that 

f Nf -Q\^dx<cf dist^(V/, SO(n))dx, 

'-'Up J, '-'Up J, 

where all inner products here are Euclidean. The theorem follows by using 
the fact that the Euclidean metric differs from the Riemannian metric by 
an 0{h^) term. 

It remains to prove (|4.10D . For the rest of this proof, exp denotes the 
exponential map of M while e denotes the exponential map of S. Let /p : 
Up^h be given by 

The metric on Upj^ relevant to formula ( |4.10| ) is /p%g- Let be sections of 
e* 3\f§ with 

= = (4.11) 

at G TpS. Let rj, v G TpS, and let 

g{s, t) = fp,h{[K^ + SX)]f(,,+sv))- 

Denote by Q the vector field along g{0, t) given by 

Q(0 = ^(0,0- 
To prove ( |4.1(J| ), it suffices to show that 

IQ(OP = tHkm' + \vf) + 0{t% 

Indeed, it is easy to see that Q(0) = 0. We denote the covariant derivatives 
of Q by Q', Q" etc. Using the symmetry of the connection. 



D_dg 
ds dt 



s=0 



f=0 



t=0 



s=0 



= -(^(0) + sx(0) + 7] + sv) 
= x(0) + v. 



s=0 
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Let R denote the curvature of g. Using the symmetry of the connection and 
assumption ( [4.11D , 



s=0 



f=0 



= k(|(0,0),Q(0))|(0,0)+|§| 



f=0 



s=0 



= + £v,,+sv(^ + sx) 



0. 



s=0 



So, 



<Q,Q)'(0) = 2(Q',Q)(0) = 0, 

(Q, Q)"(0) = 2{Q', Q')(0) + 2(Q", Q)(0) = 2 |x(0) + v\\ 
(Q, Q)"'(0) = 6<Q', Q")(0) + 2{Q'", Q)(0) = 0, 

and dODI ) follows. ■ 

This local rigidity theorem is the basis for proving the following "global" 
rigidity theorem: 

THeorem 4.2 There exists a constant C > such that for every h G (0,^o) i^^d 
every fh e W^'^(Q;,; R") there exists a section e W^'^(S; T*M|s (8) R"), suc/z that 



Jn,, 



Id//, on- 7T*q,,p rfuolg < C/z^ {£,,[/,,] + 1} . 



and 



<C {£,[/] + !} 



(4.12) 



(4.13) 



Proof: Let p;, G L^(§;T*M|s <8) R") be a section satisfying the assertion of 



Theorem ^]T|. We define 

(\h = to* [Kh Tt* (e*rf/;, o e*n o TiW (8) e*dvolg)) , 

Consider now the integral 



(4.14) 



I = I \dfh oil - n*(\h\^ rfvolg. 
Jn,, 
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Then. 



J = I 71* {\dfi, on- ncihf rfvolg) 



*@ ^71* (id/,, on- nqhf rfvolg) ((K,, o O) cD*7T*rfvolg) 
9 r {(X;, o O) d)*7T*(ivolg A e* 7i* (id/,, on- nqh\^ dvolg)) 
©iEl^ r {(^^ o O) a)*7i,dvolg A n^-e (|d/„ o n - 7i*q,p dvolg)) 

^(7i*dvolg) CO* [{Kh o O) TT* ([e*d/,, oi*n - TT* e* q^p e*dvolg)) , 

where here and in the remaining of the proof we write above relation signs 
the equation number from which the relation follows. 

Using next the fact that parallel transport is norm-preserving, adding and 
subtracting n*(D*pjj, and using the Cauchy-Schwarz inequality, we get 

I<2h+ 2I2, 

where 

h = I (7i*dvolg) (0* [{Kh o O) 71* {\edfh oi*n o n*W - n* (D*Ph\^'edv ol^)] 
^ s 

h = I (7i*dvolg) CD* [{Ku o O) TT* {\n*(D*Ph - n* e q,, o 7i*^p e*rfvolg)) . 
Consider Ji: 



h < 



(4.12) 



^ ^(7i*rfvolg) CD* {/,, TT* (fe*rf/,, oi*n o n*W - 7i*cD*p,,p e*(lvolg)} 
^ ^(7i*rfvolg) CD* {/,, 7t* (dist^(d/,„ SO(n)) oi'i^rfvolg))) 
J (7i*rfvolg) CD* |/,, 7i*e*(ivolg} 



+ 



h" 

W C ^(7i*rfvolg) CD* {^2,, TT* (dist^(rf/fo, SO(n)) oe(e rfvolg))| 

+ Ch^ J (7I*rfvolg) CD* {]<:2fe 7T*e*dvolg} = ha + hb 
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We then basically revert the steps we did before: 



Ii„ P C r (Ti^dvolg) CD* \^K2h 7i*e* (dist^(d/;„ SO{n)) dvolg)} 
C ^(Ti^dvolg) CD* {^2/, e* 71* (dist^(rf/;„ SO(n)) rfvolg)} 
^ J" CD* |]<C2fo e* 71* ^dist^(rf/;„ SO(n)) dvolg) A cD* 7T*(^volg| 
C I e* |K2;i e* 7i* (dist^{d fh, SO{n)) dvolg) A cD* 7i*dvolg| 



(0 ^ J" 71* ^dist^(rf/fo, SO{n)) rfvolg) e* {K2;, cD* 7T*dvolg} 
§ C r7i* (dist'(d/,, SO(n)) dvolg) Ch^m £h[fhl 

\J s 



and 



'lb 



= Ch^ 7i*rfvolg CD* ^Kxh 7i*e*rfvolg| 

J 7i*tivolg CD* {K21, e 7T*rfvolg} 



^ Ch- 



I- 



dvolg = Ch^ \Q.h 



Consider next l2'. 



h ^ r (7i*rfvolg) CD* {(K, o O) |cO*p„ - e* q„ o vpp 7i*^rfvolg) . (4.15) 

^ s 
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We derive the following identity: 



eqnoW^e [KhUi. ((edfh o'Tn o n*W)'e dvol^)] o W 
9 e* (D^n^ ((K,, o n)(edfh oFn o n*W)'edvol^) o W 
cD*7i*e* ((K/j o n)(e dfi, o e*n o TiW)'e dvol^'j o \p 



(4.1c) 



= cD*7i* ((K;, o 71 o e){e^dfh o e*e*n o e*7i*W) e'e*rfvolg) o W 

^ (Dijii, {{Kh °no e){e^dfh o e*e*n o e*7i*^ o ti(D*W) e^rfvolg) 

cD*7z^ (^{Kij o'eo n){Kot'edfh o KoTeH o 'iieW o tt*cd*^) H*a*e*rfvolg) 
^ (D^Tii, {{K]i o e"o n){E*a'^dfh o H*a*e*n o 7i*H*«*^) H*a*e*dvolj5) + OQi^) 
9 0* {(K,, oe)n^S*a* [(edfi, oi*n)i*dvolg) o H*a*^) + OQi^) 
^ 0* {(K,, oe)E*an^ [(^dfh o'?n)i*rfvolg) o H*« + 0{h^) 
9 {(K,, o-e)E*a*n^ ((edfh o'eUo n*W)'e dvol,)] + 0{h^) 
= 0^H* {(Kfo o-?o S-^)a*n^ ((edfn oTU o H*\p)'rdvolg)) + 0{h^) 
@ {(K,, o-^o E-^)a*n^ [(edfn o-eUo n*W)'e dvol,)] + 0{h% 



On the other hand, using the fact that 



e* 1 = e* CD* [KhUi^ dvo\^ 

— - — cD+e [Kh e Ti^rfvolg) 
= CD* {{Kh o e)e* e* TT^dvolg) 

7bl 



- CD* ((K/, o e)H*a* e* 7T*dvolg) 
QJibl ^^^^^ H"^)a*7i*'e*rfvolg) 
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we have 

cD*p;, = cD*p/,|6* ((K;, o'i'o H"^)«*7i*e*dvolg) 

© {{K,, o'?o H"^)(a;*cD*Pft)a*7T*e*rfvolg) 
|6* (^{Kh o e"o H"^)«*7i* (7i*cD*Pft e*tivolg)) . 

Thus, 

e* Qft o - a)*p,, = 
jS* |(]<:,, oe" o H"^)«*7i* ((e*ti/fo o e*n o tt*^ - tt*cd*P;,) e*t?volg)} + 0{h^), 

and by Jensen's inequality, 

|e* quoW-(D*Pi/ < 

^ f>i,d |(]<C;, oe'o H"^ o a)7T* ([e*rf/;, o e*n o tt*^ - 7T*cD*p;,p e*rfvolg)| 
CD* CD* {(Kfo oe'o H"^ o a;)7i* ([e*rf/;, o e*n o 71*^ - 7i*cD*p;,p'e*dvolg)} . 

We then use the fact that K;, o "e" o S"^ o a < Clh^^ Jzh and the local rigidity 
theorem to obtain 

\ec\hO^}- (0*p/,|^ < -j^GfG)* [lihn^ [disi^idfh, SO{n)) oi"rrfvolg)) + OQi^). 
Substituting into ( |4l^ ): 



C 



^ J" (TT+dvolg) CD* \{Kh o O) CD* CD* |/2;,7i* (dist^(rf/;„ SO(n)) o'e e*dvolg)} 7i*e*dvolg} 
^ ^(7i*rfvolg) CD* {/2;,7T* (dist^(rf/;„ SO(n)) o e'i'dvolg)} cD* {(K/, o 0)7T*e*rfvolc 
~ ^ ^(7^*^volg) CD* {/2;,TT* (dist^(rf/fo, SO(n)) o e e*rfvolg)} 
= C ^ (7i*rfvolg) CD* {^4^71* ^dist^(rf/ft, SO{n)) o 'e e*rfvolg)| , 
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and the latter is identical to 7i„, up to the fact that has been replaced by 

Putting together the estimates for Iia, l^,, and I2 we obtain 

1 



as required. 

It remains to estimate the derivative of p;,. Writing 

(\n = (D* jX;, Ti* (^{dfh o n) e*dvolg) o vp} , 

we differentiate, 

Vq,, = (D* [dKn {e{dfh o n) e*rfvolg) o 
+ (D^ [Ku Vtt* (e{dfh o n) e^rfvolg) o 
+ (D^ [Kh TT* (e{dfh o n) e*rfvolg) o Vw) . 

The second term vanishes because rr* (e'idfh ° IT) e*dvolg^ is a top-degree 
form. The third term if 0(h) because VM^ = at the zero section of TS. 

Finally, noting that 

= d{l)Ph = CO* [dK„n^ ((7T*cD*p;0e*dvolg)), 

we get 

Vq;, = CD* ^dKh Ui, {(e dfu oe*!! o tt*^ - 7i*cD*p/,) e*rfvolg)} + 0(h). 
Using once more Jensen's inequality, along with the bound for \dKh\, 



|Vq,p < j^co^ [jnn^ (\^dfn oFn o nW - h a)>,|'i*dvolg)) + 0{h^). 

Applying the local rigidity theorem and integrating over S we recover, up 
to a 1/h^ pref actor, the same bound as above. 
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Let /;, : D.h satisfy the finite bending assumption ([2.3D, and let 

q;, : TMIs — > R" be a corresponding family of sections that by Theorem ^ 
satisfies 

f \dfhon- n q,j2 dvolg < C/z^ {£,[/,] + 1} = 0{h\ 

Jq,, 



and 



J|Vq/^^^<C{£„[/,] + l} = 0(l). 



By property ( |3.3D of the push-forward operator tt* and the Cauchy-Schwarz 
inequality, 

— ^ I dist^(q;„ SO(n)) Ti^rfvolg = f dist^(7i*q;„ SO(n)) rfvolg 
l^-^/'l Js Jo,, 

< 2/1^ + 2f Id/;, o n - 7i*q,p rfvolg 

Jq,, 

= 0{h% 

from which we deduce that 



5 Compactness 



The results of the previous section can be summarized as follows: let 
G W^'^(Q/,;R") be a family of mappings satisfying the finite bending 
assumption ( |Z3| ). Then there exists a family of sections q^ G W^'^(r*M|s 
R"), such that 

I |rf/foon-7i*q,,|2dvolg = 0(/z^), (5.1) 
Jq,, 

r, 7I*(ivola 



and 

|Vq„P^^ = 0(l). (5.3) 
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^'reposition 5.1 There exists a sequence (not relabeled) % that weakly converges, 
ash ^ Q, to a limit q G W^'^(S; T*M|s ® W), namely, 

In particular, q;, — > q strongly in L^{§; T*M\g <S) IR")- 

Proof: It follows from Eqs. (|5^), (^3|), and Lemma that both q;, and 
its covariant derivative are bounded in L^(S). Weak convergence follows 
from the weak-compactness of W^'^(§; T*M|s <8) R"). The weak convergence 
of q;, to q in W^'^{S; T*M\s (2) R"), the fact that weak boundedness implies 
strong boundedness, and the Sobolev embedding theorem imply that q;, 
strongly converges to q in L^(S; T*M|s <8) R")- B 
A notational convention: we will henceforth write 

(\t = '^h° and q-^ = q o i^. 

They belong to W^'^{§; XS* ® R"). We write 

q" = q/! ° f" and q" = q o i". 
They belong to W^'\§; TS ® R"). 
CoroCCary 5.1 

lim f Irffft o n - 7i*qp rfi;oL = 0. (5.4) 

Proof: Proposition ^]T] together with Eq. ( |5.1D gives the desired result. ■ 
'Pro-position 5.2 q G SO{n) a.e. 

Proof: By the Cauchy-Schwarz inequality and the invariance of SO(n) 
under parallel transport, 

I dist^(7i*q, SO(n)) dvol,, < Ih^ lh\fh\ + 2f o n - n\f rfvolg. 

Both terms on the right hand side tend to zero as — > 0, hence 

lim f dist^(7i*q, SO{n)) dvolg = 0. 
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By the properties ([3.3D,(P^ of the push-forward operator Ui,, and the uni- 
form limit ( p7D , 

= lim f dist^(7i*q, SO{n)) dvolg 



lim — ^ I 71* (dist^(q, SO(n)) o n rfvol^) 

h^O \Llh\ Js ^ ' 



dist^(q, SO(n)) (lim- 



Ift^O |Q;, 



dist^(cj, SO(n)) rfvolg|g, 

which implies that q G SO(n) a.e. ■ 
CoroCtary 5.2 

\dfh?dvol^ = 0{\). (5.5) 



/ 

J a,, 



Proof: This is an immediate consequence of the finite bending assump- 
tion Q. ■ 

Define now the averaging operator ti* : L^(Q;,) — > L^(S), 

7i*((/)dvolg) 

7Z*({^) = -— j — , 

— 7i*avolg 

and consider the family of mappings F;, : S — > R" defined by 

hr = Tuifn). (5.6) 

It follows from ( |2.1|) and Lemma that 

lim f Fh rfvolgig = 0. (5.7) 
Js 

In the rest of this section we prove that strongly converges in W^'^(S; R") 
to a limit F G W^'^(S; R"), which is the reduced-limit of the sequence //;.. 



Lemma 5.1 



r \fh - Fh o 7ipiuo/g < Ch^j- \dfhfdvol(^ 
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Proof: By the definition of F/,, 

n*{fh - Fh o 7i) = 0. 
So, we apply the fibered Poincare inequality, 

and integrate over S with respect to the measure 7i*rfvolg. 
Troposition 5.3 

''\dF„-q^l\'dvoky=0{h'). (5.8) 



Proof: Using Corollary Lemma ^3] and the fiberwise integrability of fh, 

r ff\ 7Z^(/,,7I*rfvol3|g A CO) n^ifhco) 1 .... A r^^hW 

— 7i*(7T*avolg|g A a;) u-^co v^rf 

So, 

dFh = [n^idfh Aco) + n^{fh A dco) + (n\dniX (//'^)] ("^^d 0{h)). 

(5.9) 

By Lemma we have 

(H|5Qj,(/„a;) = 0. (5.10) 
Similarly, using also the fact that ti+o; is constant, 

ni,{da)) = d{n^(jo) + {n\da,)^ cu = 0. 

So, 

7i^(/;, A dco) = 71* ((/;, - Fh o n)dco) , 

and 

dFh - '^I = ^ [[^h - ^*^]) A a;) + 71* ((/„ - f , o 7i)da;)) (mod OQi)). 

Write dco = a Aco + where « is a 1-form and 7i*j6 = 0. Choose C such that 
|«p < C. By (|3^), 71 is a Riemannian submersion with respect to g. So, by 
Corollary we have 

l^^'^-^'Il = ^ ^ ?f ^ d^/'^ - ^) + (1//^ - ° ^1^)) (^od 0(^)). 
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By Jensen's inequality and the Cauchy-Schwarz inequality, 

m-qll' < {n. [\dfH - n*qlf a;) + Cn, (\f„ - F, o npo;)} (mod 0{h^)). 



By Lemma and equations ( [3.14D , we have 



7T*q" o (cr © t) = (cj e i)* o p o Tiql = 7i*q". 

So, by Lemma 

|rf/,-7i*p;f = |d/„on-7i*q;on|' 

< 2 o n - 7i*q,, o (cT © L)f + 0{h^) 
= 2\dfh on- 7T*q,,p + 0(/|2). 

So, integrating and using Lemma and Corollary we obtain 



= (3 + 0(h))J^ ^2 \dff, on- 7z*q,,f + CI/,, - f o Tip j rfvolg (mod 0{h^)). 

The first term on the right hand side is 0{h^) by (|5.1D. The second term is 
0(/z2) by Lemma ^ and (^). 



CoroCCary 5.3 



lim -f- IrfF/, - q"| rfuoZgig = 0. 
''^o Js 



Proof: The corollary follows from Proposition and the strong conver- 
gence of q^ to q. Proposition ^]T[ ■ 

Le-mma 5.2 Fh is uniformly bounded in W^'^(S, R"). 

Proof: By Corollary p3| , rfF;, is uniformly bounded in L^(S;R"). So, by the 
Poincare inequality and Eq. ([5.7|), F;, is also uniformly bounded in L^(S; ]R"). 
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l^heore-m 5.1 The sequence F^has a subsequence that converges to F G W^'^(S;]R"), 

Fh^F mW^'2(S;R"). 

Moreover, dF = q". 

Proof: By Lemma and the Sobolev embedding theorem, after passing 
to a subsequence, f converges in L^(S; R") to a limit f . By Corollary 
dF;, converges strongly in L^(S; R") to q". So, we conclude that f converges 
strongly in PV^'^(S;R). Since the limit must still be F, we conclude that 
dF = q". Finally, since q" G W^'^S; T*S R"), it follows that F G W^'^{§; R). 
■ 

CoroCCary 5.4 



Proof: The result follows from Lemma pTTl , Corollary Corollary |5.1| and 



To conclude, we have shown that has a subsequence that reduced- 
converges to a limit (F, q-^). 

6 Recovery sequence 



In this section we show that for every pair (F, e ^ there exists a 
so-called recovery sequence fh G W^'^(Q;,;R") that reduced-converges to 
(F, q-^), such that 




Theorem O . 



Let 



^ = W^'^(S;R") X W^'2(S;XS* ® R"). 



lim = £iim[F,q-^], 



where is given by 
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Let (F, q-^) G JT. We construct a recovery sequence /;, G W^'^(Q;,;R") as 
follows, 

fh=Fon + nq"- o A, (6.1) 
or in explicit form, for £, G 3\fS, 

= F(7i(5)) + q^^^^i^). 



^Proposition 6.1 The derivative of the recovery sequence satisfies 

dfhO o = TidF + 7i*(Vq"^) o A 
dfh o I = 71* q-^. 



(6.2) 



or in explicit form, for ^ G J^SandX = X'^eX^ G {n*T§®n*J<S)^ ^ (TSeXS)„(£), 
[dfn o (a © 0]5(X" © X^) = rff „(i)(X") + qi(£)(X^) + (Vxi,q^)„(£)(^). 

Proo/: Let ^ G and r] G (K§)„(f). Recall that 



d 



hence. 



{dfhOL)i.{ri) = -fj^{^ + tri) 
d 



t=o 



or in compact notation. 



-[F(7i(^)) + q;^(,)(5 + fr])] 

rf/;, o i = 7I*q"^. 



f=0 



Let then a:l ^ $, such that X = «(0) G (TS)„({). Recall that ^{(X) = y{0), 
where y : I ^ 3\fS is any parallel normal field along a. Then 



{dfn o c7)5(X) = - [f («(0) + qio(^(^))' 

= „(i)(X) + (Vxq^)„(,)(5) + q;^(,) [ ^ 
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The last term on the right hand side vanishes because y is parallel. So, in 
compact notation, 

dfi, 00 = TidF + 7T*(Vq"^) o A. 



(Proposition 6.2 fu reduced-converges to (f, p-^), namely, 



lim f I f/j - f o Tip dvoL = 



lim f \dfh on- nidF © q^)|^ dvol^, = 0. 



Proof: For every ^ G Q;„ 

ifn - F o nm = = OQi), 

hence 

f \fh - F o dvolj = OQi^). 
Jn,, 

Next, by the Cauchy-Schwarz inequality, 

\dfh on- n{dF © q^)p < 2 o (a © i) - n*{dF © q^)p + 2 o (n - a © i)P 

< 2|7i*(Vq^) o Ap + 2\dfhf |n - a © tp. 

Since by Lemma ^Aj \Tl - o ® i\ - OQi) and dfh is uniformly bounded, i.e., 
satisfies (|575|), it follows that 



f MM 

Jq,, 



|n - cT © dvolg = 0(/z2). 



Finally, since uniformly for every ^ G Q;,, 

|H (Vq^) o AU) = \{Vq^Um = 0{h), 

it follows that 



f |7i*(Vq^)o Aprfvola = OQi^). 

Jq,, 
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Lemma 6.1 ¥ or all A G GL{n), 



dist(Id +A, SO(n)) 



< C min(|AUA|2). 



Proof: The 0(|A|) bound follows from the fact that all the terms on the left 
hand side are 0{\A\). The 0(|Ap) bound follows form the fact that |A + A^|/2 
is the first-order Taylor expansion of dist(J + A, SO{n)) at A = 0. ■ 



Proposition 6.3 



lim £;,[/;,] = £/,-,„[f,q-^]. 

h->0 



Proof: Note first that by Lemma ^]T] and Proposition |63t 

dfhoTl = dfh o (cj © - ^^//, o cj o 71* II o(A ® TiT") + 0(/i2) 
= 71* (df © q^) + 7i*(Vq^ -dPo II) o (A 7i*P") 
+ (|Vqi + |d/,|)0(/z2) 

By the invariance of SO{n) under parallel transport, 

dist(rf/;„ SO{n)) = dist{n{dF © q^) + 7z*(Vq^ - dF o II) o (A nP"), SO(n)) 
+ (|Vq^| + |d/„|)0(/z^). 

Hence, by the uniform L^-boundedness of d/;„ 

lim £,,[/fo] = lim -^i- dist^(7i*(df ©q^)+7i*(Vq^-df oII)o(A(87zT'I), SO(n)) dvoL 

Note that for ^ G Q;„ 7T*(Vq^ - if o II) o (A tz^P")!^- = 0(/z), hence 

lim h^^i,[fh]= lim I dist{n{dF ® q^),SO{n))dvoL 

=£dist{dF © q^, SO(n)) dvol^i^, 

and therefore if df © q-^ ^ SO{n) a.e., then 

lim £;,[/ft] = oo = £iim[f,q-^]. 
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It remains to consider the case q = dF © q-*- G SO{n). We have 

dist(7z*q + 7T*(Vq^ - dF o II) o (A P"), SO{n)) 
= dist(Id +71* (Vq^ - dF o II) o (A (8 P") o nq-\ SO{n)). 

It follows from Lemma |]T] that for A e GL{n) with \A\ = OQi), 



dist^Id +A,SO(n)) 



Take 



A = 7z*(Vq^ - df o II) o (A 7i*P'l) o 7Z*q"^ 
or more explicitly, for v eM" 

Using the fact that A^ = O(^), we obtain 



(6.3) 



lim £;,[/fe] = lim -^i- 

= lim ^rr^t 



A+A^ 



dvoL 



tr(A^A + A2)dvolg. 



(6.4) 



At this point, it is helpful to introduce index notation to clarify the sense 
in which tensors with several upper and lower indices are composed and 
transposed. Let p e S, and let Xi, . . . , X„, be a basis for TpM such that 
Xi, ...,X„, is a basis for TpS, and X„,+i,. . .,X,„+„, is a basis for KSp. To 
keep our notation concise, we use the following convention for ranges of 
summation: 



a, b, c, 

u, V, lU, 



run from 1 to n, 
run from 1 to m, 
run from m + lton. 



The orthogonality of TpS and KSp implies that Qau = 0. 

Let Yi, . . . , y,i, be a basis of ]R". We reserve Greek letters for indices associ- 
ated to the Y's. For ^ G XSp and 7] G R", write 
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So, for example, 

A,rj = (V„q«-q«l4)(q-i);rr]^y«. 
By equation and Lemma we get that 

lim Ehifhl = f W rfvolgig, 

where 

W = lim f-l-7i*(tr(A^A + A^))] . 
We now calculate W explicitly. Define the section M of 3\fS (2) 3\fS by 



1 

M = lim -tTcJA (2) A). 



At the point p e S, we have 



+ -{VX - q?l4)(q-0;(Vfcq? - q^, II^,)(q-0«M 

Since we have chosen Q/, to be the set of all points of distance less than h 
from S, we have M"" = kq"^, where k is the volume of the k-1 dimensional 
unit sphere divided by + 2. Hence the reduced energy density is 



+ ^9"^((q-')aV„q2 - lt)((q-^);V,qf - 



(6.5) 



Since (dF)^ = q", it follows that 

= V«q«. (6.6) 
Differentiating the equation = P" o = pH o q"^ o q o i-"- = p'^ o q-"-, we obtain 

= VP'I o q^ + Pj o Vq^, 
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or, in index notation, 

(V.p;:)(q-^)^ = -(V«(q-i)^)q;;. (6.7) 

By equations and (|6.7|), the fact the metrics are covariant constants 
and the orthogonality of q, 

Wu){c^-'t = -C^aic^-XK = -e«,S3'^(V„qr)q« 
= -e«,g'"^(V,qr)q« = e,,g''^(V,q^)qr. 

Substituting this last identity as well as the symmetry of the second fun- 
damental form into the second term of (|6.5D, 

9"*((q-')^V,q« - lt)((q-i);V,qf - 1^) = 

= 8"^(ea,8^^(V,q«)qr - g''^g„jt)((q-i);V,qf - If,,) 
= 9'"e„,g^^((V,q«) - q«lt)(qr(q-i);V,qf - qllf,,) 
= |P"o(Vq^-qlloII)|2. 



Thus, 



W = J|Vq^ - q" o II |2 + lipj o (Vq^ - q" o II)|2. 



2' " " '2' 
Noting that I ■ P = IPj o -I + o -I, pj o q" = Id, and P^ o q^ = 0, we have 

W = k\PI o Vq^ - Hp + ^\P^ o Vq^p. 

The proposition follows immediately. I 



7 Lower semicontinuity 

In this section we show that a sequence /;, that satisfies the finite bending 
assumption ( [2.3[ ) and whose reduced-limit is (F, q-"-), satisfies the lower- 
semicontinuity property, 

liminf £,,[/,,] > £iim[f,q^]- 

We first pass to a subsequence /;,^ so that 

lim ^h,[fh] = liminf £;,[/;,]. 
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Thus, in the following arguments, we may freely pass to a further subse- 
quence; from now on, we drop the subscript k. 

Let qh be a sequence obtained from the rigidity Theorem and let q/j 
be an orthogonal projection of q;, on SO{n), i.e., q;, G SO{n) and Iq^, - 
qh\ = dist(q;;, SO(n)). Applying the Cauchy-Schwarz inequality and the 
invariance of SO{n) under parallel transport, 

\nqh - nc\hf < 2 \dfn oU- nquf + 2 dist^(rf/;„ SO (n)). 

Averaging over Q;„ using once more the properties of the push-forward 
operator ti*, estimate (pTTD, the finite bending assumption (|2.3D, and the 
estimate of the volume form discrepancy. Lemma we have 

f\q,-~qn\^dvo\,i^=0{h^). (7.1) 
Js 

It follows from Eqs. and (^ that 

\dfh on- 7i*q,,|2 dvolg = 0(/|2). (7.2) 
Passing to a subsequence and using the L^-convergence of q^ to q (Propo- 



sition 5.1), we have 



lim f |q;, - qNvolcis = 0. 

Note that we can only guarantee the convergence of q^ to q in L^(§; r*M|s (8> 
R"), whereas q^ converges to q weakly in W^'^{§; r*M|s <8) IR"). The reason 
for defining the possibly less regular sequence q^ will be made clear further 
below. 



^Proposition 7.1 The sequence 



% = — — {7.3) 



has a subsequence (not relabeled) that weakly converges, ash ^ 0, in L^(S; T*M|s<8) 
R"); ive denote the limit by a. 

Proof: This is an immediate consequence of Eq. ([7.1|), as it follows that the 
sequence is bounded in L^(S; T*M|s <8) IR")- ■ 
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^Proposition 7.2 Let F/, = TLi,{f}i) define 

fh=Fhon + nq^ o A. 

Then 

lim f \dfh - dfh\^ dvolc^ = 0{h^). 



Proof: Compare the definition of fij with that of the recovery sequence ( |5TTD . 
By an argument similar to that used in Proposition we find 



dfhoo = ndFh + 7i*(Vq,f ) o A 
dfh o L = n(\^. 

Using the invariance of the inner-product under parallel transport, along 
with the Cauchy-Schwarz inequality (twice) and Lemma pTT] , 

\dfh-d0 = \dfn on- dfh oUf 

< 2 \dfh on- dfya e + 2 \d0\{a ® l) - Hp 

< 4 \dfu o n - n{dFn ® q,t)P + 4 |7z*(Vq,^) o Ap + 2 \dfif OQi^). 

We then average over Q;,. The first term on the right hand side is 0{h^) 
by Eq. (pTTl), Proposition |53| , and Lemma The second term is 0{h^) 
because Vq^ is bounded in by Eq. (|5!3|) and |A| = 0{h). Finally, the third 
term is 0{h^) because is uniformly bounded in L^(Qfo) as obtained by 
combining Eq. (p^). Lemma |XT| , Lemma Eq. ( |5^ ) and Eq. ■ 

Troposition 7.3 Let &e defined as above. Then the sequence 

dfh - dfh 

hh = (7.5) 



weakly converges to zero in the sense of Definition p7]\ . 

Proof: Let O G 0^(0;,^; TD.h, Denote by d* the dual of the exterior 

derivative with respect to the pairing of g. By Lemma d* is also 
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the dual of the exterior derivative with respect to the pairing of /i*g. 
Integrating by parts. 



Using the Cauchy-Schwarz inequality and Lemma 

2 , 

I iKQ)Wdfh-dfy,0)dYohQ <C'(0) I |/,,-/^po^;,dvol. 



<C(0)|£ |/,-/„|^rfVOlgj. 

Using the fact that ni,{fh) = ni,{fh) and applying the fibered Poincare in- 
equality, 

C(0) 



< 



< 



C(0)/z2 



P/jI Js 

,2 



71* (l/Zi -/^P) 7I*dvolg 

n^{\dfh - d0) Ti^dvolg 



= C(0)/z^ f \dfh - dfh\^ dvolg. 



Dividing by and applying Proposition \772\ , 



lim 

h->0 



8) 



dfh- dfi 
h 



dvolf 



< C(0) lim f It?/,, - rf/,,r dvolg 
= 0. 



Since C°°{Clh,; Tdh^^W) is dense in L\Oh^; F (81R"), this equation holds 



for all O G L2(Q. ; rQ;,„ ® R"). 



'Proposition 7.4 The sequence 



h 



(7.6) 



/zfls subsequence (not relabeled) that weakly converges, ash ^ in L^(S; r*S (8> 
R"); we denote the limit by c. 
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Proof: This is an immediate consequence of Proposition [5.3|, which together 
with Lemma implies that 



\dFn-qlrdvoki,=0{h% 

s 

hence is a bounded sequence in L^(§; T*S (2) R"). ■ 

We now turn to estimate Eh[fh] as h — > 0. By the invariance of SO{n) under 
parallel transport and the fact that G SO{n), 

dist^idfh, SO{n)) = dist^(d/,, o n, SO(n)) 

= dist^(d/;, o n o n%\ SO{n)) 
= dist^ {I + h Gh, SO{n)) , 

where G;, G L^(Q;,; GL(n)) is given by 

dfhoTlo 7i*q"^ - Id 
Gn = — . (7.7) 

Thus, 

Wh] = dist'(Id+/zG,„SO(n))rfvol3. 

We start by making a few observations about G^: 

'Proposition 7.5 The sequence Gh g L^(Q;,;]R" ® R") zs hounded, namely, 

f \Gh\^dvol^ = 0(1). 
Jq,, 



Proof: This is an immediate consequence of the estimate ( [7.2D , as 



f |G;,P rfvolg = Id/;, o n o n%'^ - Id p dvolg 
= ^£ o n - nqhfdvol, = 0(1). 
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(Proposition 7.6 Gh can he expressed in the following form: 

Gh = \oUo 7i*q-^ + n\{ah + c,, o p") o q-^) 

t^/^ o (n - cj e ( + a o 71* II o(A 7i*Pll)) _i 
+ 1 ° (^*% ) 

+ 7i*(Vq,-^ - dFu o II) TiT" j o Tf%^ 

- 7i*(Vq^) o |a (2) 71* II o 1^ (2) 7T*P" jj o 7i*q;;"\ 

Proof: This follows after straightforward algebraic manipulations from the 
definitions ( p^ ) and ( [73D of a;, and b/,, formula for rf/;, o (cr © i), and the 
definition (Ob of C;,. ■ 



(Proposition 7.7 Let 

im - hi 



hVJh\ = — t Xh dist (Id+/iG;„SO(n))dz;oZg, 

Q,, 



where 



and let 



Then, 



Xh 



1 IGfol < 
1 otherwise, 



Jhifh] =j'Xh 



Gh + GT 



dvoky 



lim(/,[/,] -£,[/,]) = 0. 



Proof: By definition. 



Jhifh] - ^h[fh] = f Xh (al - jS^) rfvolg. 



where 



dist(Id +h Gh, SO{n)) 
ah = and = 



Gh + Gl 
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Using Lemma ^]T] with A replaced by h. G/„ we obtain 

\ah-^h\ < Cmm{\Gh\,h\Gh\^). 

When the indicator function Xh is non-zero, \Gh\ < h~^^^, hence 

\al - < \an - j6;,| {\an - + < Ch\GH\^ ■ 2\Gh\ = 0{h^' 
from which follows that 

MM - h[fk]\ < f Xh\al - Pi\ dvo\ = OQi^l\ 
Jn,, 

which completes the proof. 

Since, trivially, £;,[/;,] < 8,h[fh], if follows that 

liminf Eh[fh] > liminf /;,[/;,], 

hence it only remains to show that 

liminf /,,[/,,] > £iun[f,q^]- 

To this end we write Gh = G|^^^ + G^^\ where 

Gi'^ = G„-^-7il(copll + a)oq-i], 



and 



Gf ' = ^ + 7il(c o P" + a) o q-i], 
" h 



where A is given by Eq. (|63|), namely, 

A = 7i*(Vq^ - q" o II) o (A 7t*P") o 7z*q"\ 

Thus, 

h[fk] = \( Xh\G';^ + Gffdvok 
^ Jq,, 

+ ^X.e(Gi^^ + G^)^Gf + Gf'),dvol, 
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and 



limmf /,[/;,] > limmf i£ Xu \Gf + Gf\^dYo\ 

+ lim inf \£xk^ [g^;^ + Gf^\ Cf + ^) rfvol^. 



Troposition 7.8 



lim inf \ f Xh |GP + Gf^^l'^dvoL = lim inf I f |GP + Gf^^\^dvoL. 

Proof: Consider the difference 

A„ = f (l-X„)|Gf + Gf>dvolg. 

Jq,, 

Using the Cauchy-Schwarz inequality (twice), the uniform bound \ A/h\ < 1, 
and Lemma 

A,,<4f (l-x;,)|Gfpdvolg 

< 8 f (1 - x,)7z* (iVq^^ + „ pii + ^)|2 + „ jj)|2\ ^^^1^ 

Jq,, 

< C r ((1 - X/,) o ^n) 71 (iVq^p + |(c o P" + a)p + |(q" o II)|2) a; A 7]. 

(7.10) 

Using Lemma and Proposition 

I ((1 - X/;) ° f ;<)t^ A 7] < C'l (1 - Xh)dvo\o^ 

< Ch^'^f |G,/dvolg = 0(^1/2). 
Jq,, 

It follows that (1 - Xh) ° fi/i is the indicator function of a set of measure 
tending to zero with h. So, since n* (|Vq-^p + |(c o pH + a)p + |(q" o II)p^ is 
integrable in Q^g, inequality ( [7.1UD implies that tends to zero as ^ — > 0, 
which completes the proof. ■ 
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'Proposition 7.9 



liminf if \Gf + GfUvol,>Eum{'P,(\% 



Proof: The integral 



depends on the sections a and c. It is easy to see that for every finite h> 
this integral is minimal for a, c = 0{h^), hence 



and the last identity was proved in Section ^. 



liminf 



dvoL = SiimK q-^]. 



Proposition 7.10 



liminf 



= lim inf 



9 



G^ + Gf,Gf^ + Gf)dvok. 



V. = r (1 - X/Oe ( + Gf, Gf + Gf ^) rfvolg. 



Proof: Consider the difference 

A,, 

Using the bilinearity of g and the Cauchy-Schwarz inequality, 

|A,|2<4^ |G;;'NvolgJ|£ \XhGf^\^dvol,^. 

The first term on the right hand side is uniformly bounded by Proposi- 
tion whereas the second term tends to zero by the same argument as in 
the proof of Proposition ^!8| . ■ 
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'Proposition 7.11 



lim - f efcf) + Gf,Gf^ + cf] dvol, = 0. 
Proof: Note that 

e (g1^) + Gf, Gf + Gf = g [{G^ + gJ^)') o n^,, (Gf + Gf o n^,) . 



Using Proposition the sections G|^^^ o 7i*ci can be rewritten as follows: 

G^l^ o 7i*q = bfo o n o n{(\~^ o q) 

+ n[{ah + Ch° f ") o Q;^^ o q] - 7i*[a + c o p"] 

+ 7i*Vq,-J- ° 1^ 7i*P" j o n{c\l^ o q) - 7z*Vq-^ ° *^ nP^^ 

- n{dFh o II) 7i*P" j o 7T*(q;;^ o q) + 7i*(q'l o II) « ® tiT" j 

rf/^ o (n - cr e i + cr o 7Z* II o(A ® tiT")) _^ 
+ ° ° Q) 

- 7i*(Vq^) o |a (8) 71* II o 1^ (2) 7i*P" jj o n{c\~'^ o q). 

Consider the first line of the right hand side: b;, weakly converges to zero 
in the sense of Definition (Proposition whereas is bounded in 
L°°{§;W (8 TMIs) and strongly converges to q"^ in L^{§;W ® TM\s). It is a 
known fact that the product of an weakly converging sequence and an 
strongly convergent sequence that is bounded in L°° weakly converges 
in to the product of the limits. We therefore conclude that the first 
line weakly converges to zero in the sense of Definition By a similar 
argument the second, third, and fourth lines weakly converge to zero as 
well. For the fourth line, we use also Proposition The fifth line is 0(h) 
by Lemma whereas the sixth line is 0(h) due to the A-factor. Thus, 
G^^^ o 7T*q weakly converges to zero in the sense of Definition PTT| . Similarly, 
Gj^ o 7i*q weakly converges to zero. Noting that 

(Gf + GfVn*q = 71*^1 + 71*^2°^ 



for suitable ^1,^2, we apply Lemma ^ to obtain the desired result. 
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CoroCtary 7.1 



liminf ]u{fh\ > E,iu„[F, (\^]- 



Proof: This is an immediate consequence of Eq. (^!9|) and Propositions 



CoroUary 7.2 



liminf £;,[/;,] > £;,>„[f, q-^]. 



8 Examples 

8.1 Plates and shells 

Plates and shells correspond to the case of n = 3 and k = 1, that is, the limit- 
ing manifold § is a two-dimensional surface. For k = 1, the limiting energy 
functional further simplifies, as the second term in ( |2.6| ) vanishes. This 
is because for k = 1, 

o Vq^ = 0. (8.1) 
Indeed, let u be a local unit length section of KS. Then, Vu = because 

= dQ{u, u) = 2q{u, Vm) 
and ?\fS is one-dimensional. So, 

(Vq^)(t/) = d{c\^{u)) - q^iVu) = d{q^{u)). 
On the other hand, since P^ is an orthogonal projection, and q is orthogonal, 
g(P;^ o (Vq^)(i/),w) = e((Vq^)(i/),q^(i/)) = 

= e(c^(q^(M)),q^(w)) = ide(q^(u), q^(u)) = 0. 

Since u spans 'NS, equation { ^7\^ follows. 

Finally, we note that since q G S0(3), it follows that q-"" is unambiguously 
determined by dF, which means that the f : S — > R'^ fully characterizes 
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1 /3, the limiting functional for 

IlNvolgi,. (8.2) 

This integral has a well-defined physical meaning. It is the mean square 
difference between the induced second fundamental form PjJ o Vq-"- (the 
gradient of the normal to the surface) and the intrinsic second fundamental 
form, i.e., it is a bending energy. 

The limiting energy ( p^ applies equally to plates and shells, the only differ- 
ence being whether II = (plates) or II (shells). It also applies equally 
to Euclidean and non-Euclidean plates/shells. The plate/shell is Euclidean 
if there exists an immersion F : S — > R'^, such that F*e = g|s andPlloVq-"- = II, 
or equivalently, if g|s and II satisfy the Gauss-Codazzi-Matnardi equations 
[ po 76[ j. If Qs and II are incompatible, then the plates/shell is residually- 
stressed and the mtnimizer of ( p.2D has non-zero energy. 



the limiting configuration. Noting that k = 
plates and shells is: 

£iim[F] = ^J|P!oVq^- 



8.2 Rods 

Rods correspond to the case n = 3 and k = 2, which means that S is a one- 
dimensional submanifold. For n - A: = 1 it is dF that is uniquely determined 
by q^. 

The remarkable fact in the case of a one-dimensional submanifold is that 
there exists a limiting configuration, (F, q-"-), such that the limiting energy 
Slim vanishes. That is, dF ® q-^ e SO{n) and 

o Vq^ = and Pj o Vq^ = II . (8.3) 

This means that there exists a sequence of approximate minimizers for 
which 8,h[fh] = o(l). Hence, to obtain a finer limiting structure one should 
divide the physical elastic energy by a higher power of h. 

We now show that a solution satisfying ( ^3] ) does exist. In particular, we 
show that equations (|8.3D are equivalent to a system of linear ordinary 



differential equations. Let p e §. Since q" o P'^ = Id -q-^ o P^, it follows that 
q" o Pj o Vq^ = Vq^ - q^ o P;^ o Vq^. 



56 



Substituting using equations ([5^), we obtain 

Vq^ = q"oII. (8.4) 

On the other hand, equation dST^) immediately implies equations ( |8^) . 

Moreover, an equation for ql' can be derived from ( |8.4| ) as follows. Since 
for every ^ G 3\fS and rj eT§ 

e(q^(aq"(r])) = 0, 
differentiating and using equation (|83^), we obtain 

e(q^(aV,,q"(r])) = -e(V„,q^(aq"('?)) 

= -e(qll o II(rj', ^), qll(r?)) = -g(II(7]', e), r]). 

Introducing the metric adjoint of the second fundamental form II : TS X 
T§ — > KS, we have 

e(q^(^),V,,qll(7])) = -g(<^,II(T],Tj')) = -e(q^(aq^ o n(rj,7]')), 
from which we conclude that 

Vq" = -q^ o II. (8.5) 

Since S is one-dimensional, (|8.4D and ( ^3| ) form a linear ordinary differ- 
ential system, which can be solved globally. In fact, they are generalized 
Serret-Frenet equations, which uniquely determine the shape of the limit- 
ing curve. 



9 Discussion 

The main contribution of this paper is to unify the analyses of dimensionally- 
reduced elasticity models, resulting in a limiting model that covers a va- 
riety of known cases as well as cases that have not yet been rigorously 
treated, such as non-Euclidean shells and rods. The generalization has 
been attained by performing the entire analysis within the framework of 
Riemannian geometry. While being at times technically cumbersome, the 
Riemannian formalism is the appropriate framework when considering 
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incompatible elastic bodies, so that the notion of a Euclidean reference 
configuration becomes irrelevant. 

The entire analysis rests on the so-called finite bending assumption, whereby 
there exists a family of mappings /;„ such that ^hifhl = 0(1). As noted by 
Lewicka and Pakzad [ p!.P10[ j, this condition is equivalent to the existence 



of a W^'^ isometric immersion of (S, g|s) in R". The equivalence between 
the two conditions in an immediate corollary of our results: if the finite 
bending assumption holds then we have shown that has a subsequence 
that reduced-converges to a W^^(S;R") isometric immersion. Conversely, 
if a W^'^ isometric immersion exists, then the recovery sequence defined in 
Section ^ satisfies the finite bending assumption. 

This work raises a number of questions of interest. 



Under what conditions on Q\g does the finite bending assumption 
hold? There exists a large amount of work on Holder regular or 
smooth isometric immersions in Euclidean space. However, we are 
not aware of similar work for Sobolev maps. 



The proof of local rigidity. Theorem ^^11 , relies on the rigidity theorem 



proved in [ j^JMU2| ] for mappings R" — > R". Generalizations to map- 



pings between Riemannian manifolds would be of much interest. 



How does ff, approach the limit (f, q-^), and in particular, how does 
£;j[/ft]~£iim[f/ Q"^] scale with/z. The (non-rigorous) analysis in [ ]ESK09a[ J 



indicates that the deviation of fh from the limit F is focused in a bound- 
ary layer of width 0{h^^^) in the vicinity of the boundary of S, and 
consequently, E,h[fh] - £iim[f,q"^] = 0{h^'^). 
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